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Abstract 

We present the results of a lattice QCD calculation of the form factor relevant 
to B -» K*j decay. Our final value, T(0) = 0.24 ± 0.03±g;gf, is obtained in the 

quenched approximation, and by extrapolating M^ 2 x T H ^ K * (q 2 = 0) from the 
directly accessed i?-heavy mesons to the meson B. We also show that the extrap- 
olation from B — > K*-y* (q 2 ^ 0) to B — > K*^ (q 2 = 0), leads to a result consis- 
tent with the one quoted above. On the other hand, our results are not accurate 
enough to solve the SU (3) flavor breaking effects in the form factor and we quote 
T B ^ K *(0)/T B -*P(0) = 1.2(1), as our best estimate. 



PACS: 12.38.Gc, 13.25.Hw, 13.25.Jx, 13.30.Ce, 13.75Lb 



1 Introduction 



The flavour changing neutral decays, B — > V'j (V = K*,p,u), are induced by penguin 
diagrams. Their accurate experimental measurement gives us information about the heavy 
particle content in the loops, and thus might be a window to the physics beyond the 
Standard Model (SM). This is why a huge amount of both experimental and theoretical 
research has been invested in studying these modes over the past decade. 

The experimenters at CLEO were first to observe and measure the B — > K*^ decay 
rate [1]. Averaging over the neutral and charged £>-mesons they reported 

B(B -> K*j) = (4.5 ± 1.5 ±0.9) x KT 5 . (1) 

Today, the unprecedented statistical quality of the data collected at the 5-factories made 
it possible to do precision measurements separately for B° and B ± decays, namely 

BaBar [2] 
Belle [3] 

(2) 

BaBar [2] 
Belle [3] 

Besides, the first significant measurements of B — > p(uj) , ~f [4], opened a discussion on 
the possibility of constraining |14d/Vt s |, thus providing an alternative to the constraint 
arising from the ratio of the oscillation frequencies in the B s — B® and B d — B® systems, 
Am Bs /Am Bd [5]. 

However, when looking for non-SM effects in these decays, one should be able to con- 
front the above experimental results to the corresponding theoretical estimates within the 
SM. 1 As usual, the main obstacle is a lack of good theoretical control over the hadronic 
uncertainties. The hadronic matrix element entering the analysis of these, electromagnetic 
penguin induced, decays is 

(VV; e x )\T^(0)\B(p)) = e* a (p', A) x , 



b (b° -> ir° 7 ) 



B (B + -> ir + 7 ) 



(3.92 ±0.20 ±0.24) x 10~ 5 
(4.01 ±0.21 ±0.17) x 10~ 5 

(3.87 ±0.28 ±0.26) x 10~ 5 
(4.25 ±0.31 ±0.24) x 10~ 5 



2p° 



Ml 



m v „ \ m , 2\ , M B -m 2 v 



% W) + 



qfsT 2 (q 2 ) 



+ 



^Wa T 2 (g 2 )-T 1 (g 2 )± 



M 2 B -m 2 v 



(3) 



where q = p — p', the tensor current T^ v = isa^ u b for V = K*, and T^ v = ida^b for V = p, 
with = |[7 M , ju]. The above definition is suitable for the extraction of the form factors 
from the correlation functions computed on the lattice. The form factors Ti^^{q 2 ) are the 



1 A review on rare i?-decays, containing an extensive list of references can be found in [6]. 
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same as those computed by the QCD sum rules [7,8]. For the physical photon (q 2 = 0), 
the form factors Ti(0) = T 2 (0), while the coefficient multiplying T 3 (0) is zero. 

In this paper we show that the strategy which we previously employed to compute 
the B — > Ti semileptonic form factors [9] can be used to compute the radiative decays as 
well. Although the attainable accuracy is quite limited, we believe the values we get are 
still phenomenologically useful. In what follows we will show how we obtain T B ~* K (0) = 
0.24(3) (^i), from our quenched QCD calculations with 0(a) improved Wilson quarks at 
two lattice spacings. We also obtain T B ^ K * (0)/T B ^ p (0) = 1.2(1) although that result is 
unstable when applying different strategies and using different lattice spacings. 

2 Methods to approaching q 2 = and Mb 

Even though we work with ever smaller lattice spacings ("a"), we are still not able to work 
directly with the heavy 5-quark. Instead of simulating a meson with Mb = 5.28 GeV, we 
compute the form factors with fictitious heavy-light mesons (H) of masses Mb > Mh > 
M D) and then extrapolate them in 1/M H to l/M B) guided by the heavy quark scaling 
laws. Alternatively one can discretise the nonrelativistic QCD (NRQCD) which basically 
means the inclusion of l/(amf>)-corrections to the static limit. This, in the lattice QCD 
community, is known as the "NRQCD approach". Finally, one can build an effective theory 
that combines the above two, which is known as the "Fermilab approach". Each of the 
mentioned approaches has its advantages and drawbacks. They were all used in computing 
the B — > n semileptonic form factors and the results show a pleasant overall agreement 
(see e.g. fig. 3 in the first ref. [10]). 

Concerning the methodology employed while working with propagating heavy quarks, 
one should keep in mind that the form factors are accessed for all q 2 G [0, (Mh — my) 2 }. 
Only after extrapolating to M B) at fixed values of v-p', the g 2 -region becomes large and 
the form factors are shifted to large g 2 's. 2 The assumption underlying this extrapolation 
is that the HQET scaling laws [11] remain valid when E = v-p' > my. In the case of 
B — > 7i£u, it appears that this assumption is not particularly worrisome as the form factors, 
after extrapolating to M B [9, 12], are consistent with those obtained by the effective heavy 
quark approaches [13,14], at large g 2 's. If one is interested in the form factor at q 2 = 0, 
then in order to extrapolate from large g 2 's one has to make some physically motivated 
assumption about the g 2 -shapes of the form factors. 

Otherwise, when working with propagating heavy quarks, one can extrapolate the form 
factors directly computed on the lattice at q 2 = in 1/Mh to 1 /Mb • The useful scaling 
law relevant to this situation was first noted in the framework of the light cone QCD sum 
rules (LCSR) [15], then generalised to the large energy effective theory in ref. [16], and 
finally confirmed in the soft collinear effective theory [17]. The underlying assumption in 
this extrapolation is that the scaling law would remain valid even when the light meson 
is not very energetic (in the rest frame of the heavy, the q 2 = point corresponds to 
E = (M 2 H + m 2 v )/2M H ). 

2 v stands for the heavy quark (meson) four-velocity, so that q 2 = M| + m v — 2MbV-p'. In the rest 
frame of the heavy meson, E — v-p 1 is the energy of the light meson emerging from the decay. 
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In ref. [9] we showed that the results for the B — > nlv form factor at q 2 = 0, obtained by 
employing either of these two different ways of extrapolating to Mb, are fully compatible. 
The method of extrapolating in 1/M H at q 2 = fixed is particularly useful for B — > K*(p)^y, 
where the main goal is to compute the form factor when the photon is on-shell, T(q 2 = 0). 
This is what we do in this paper. As a cross-check of our result we also use the standard 
method (extrapolating in l/M H prior to the extrapolation in q 2 , down to q 2 = 0). 

3 Raw Lattice Results 



The form factors are extracted from the study of suitable ratios of three- and two-point 
correlation functions, namely 



Rc,tiv\ty) — 



Ci%(tx,t y ;q,p H ) 



§E C -(W') C%{t x -t y J H ) 



(4) 



a=0 



where p ' = Ph — q. The correlation functions and their asymptotic behavior are given by 



CP H (t,f) = £ e**(P 6 (z,t)Pt(0)) — ^t/ EHt 

X 



C$(t,p>) = £ e*'*(V a frt)Vi(0)) ^ 1 £^e* a {j/ t X)e p {p' t X)e 



-E v t 



4%(t x ,t y ;q,p H ) = J^^-^HVdOKMPM 



tx>t y >0 



^ lS^ (P '' A) e ~ EVtV X ( V (p'^)\T,Am(PH)) x ^ e-JM*.-*.) 



2£ 

AEyEjj 



(5) 



For the interpolating field we choose P 5 = qi'j^Q and = g7 M g, with g and Q being light 
and heavy quark respectively. We also defined \J Z H = \{0\P 5 \H)\, and \/Z~^e* jl {p l ', A) = 
| (0| A)) | . The hat over the tensor current indicates that it is 0(a) improved and 

renormalised at some scale /x, i.e., 



T^M = Z K T ] \g\, fi) l + b T (gl)am) iQa^q + ac T (g%) (d^Q^q - d^Q^q) , (6) 

where am q ^ = (1/k 9 |q} — l//t cr )/2, andm = (m q +mQ)/2. The renormahsation constant, 

Z^t 1 (do ' /•*) an< ^ the operator improvement coefficients, &r(<7o) anc ^ ^(fi'o)' are specified in 
table 1, where we also give the basic information about our lattices. By using the standard 
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Set 1 24 3 x 64, (3 = 6.2, c sw = 1-614 

200 configs; a" 1 = 2.7(1) GeV; t x = 27 
k q 0.125; 0.122; 0.119; 0.115 

K q 0.1344; 0.1349; 0.1352 

= 0.876(2); b hpt = 1.22; c T = 0.06 
Set 2 24 3 x 64, (3 = 6.2, c sw = 1-614 

200 configs; a" 1 = 2.8(1) GeV; t x = 31 
0.128; 0.125; 0.122; 0.119; 0.116 
K q 0.1344; 0.1346; 0.1348; 0.1350; 0.1352 
Zf = 0.876(2); b^ pt = 1.22; c T = 0.06 
Set 3 32 3 x 70, (3 = 6.45, c sw = 1-509 

100 configs; a" 1 = 3.8(1) GeV; t x = 34 
k q 0.1285; 0.125; 0.122; 0.119; 0.116; 0.114 
Kg 0.1349; 0.1351; 0.1352; 0.1353 

Zf = 0.883(2); = 1.20; 4 pt = 0.02 

Table 1: Details on the lattices used in this work including the values of the Wilson heavy (kq) 
and light (K q ) mass parameters; 0(a) improvement coefficient of the action csw [18] and of the 
tensor current ct [20], bx; renormalisation constant Z^\\/a) [19]. When a nonperturbative value 
is not available, we take its estimate in boosted perturbation theory ("bpt"). 

definition of the vector current form factor 

(V(p', \)\V»(0)\B(p)) = e^elip', X) Pct p' p * V{ f^ , (7) 

one can easily see that the improvement of the bare tensor current leaves the form factor 
T 2 (q 2 ) unchanged, whereas the form factor Ti(q 2 ) is modified as 

TT pr -(<? 2 ) = Uq 2 ) - ac T q " VV) ■ (8) 

Mb + Tny 
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To study the form factors' g 2 -dependence we considered the following 12 combinations 
of ph, and q: 

p H = (0,0,0) & ?e{(0,0,0) i; (l,0,0) 4 ;(l,l,0) 6 ;(l,l,l) 4 ;(2,0,0) 4 }, (9) 

p H = (1, 0, 0) & <f e { (0, 0, 0) 4 ; (0, 1, 0) 12 ; (0, 1, 1) 6 ; (1, 0, 0) 4 ; (1, 1, 0) 12 ; (1, 1, 1) 12 ; (2, 0, 0) 4 }, 

in units of (2ir/La), the elementary momentum on the lattice with periodic boundary 
conditions. The index after each parenthesis in (9) denotes the number of independent 
correlation functions Cj^ u (t x ,t y ; q,pn), for a given combination of pn and q. Those are 
deduced after applying the symmetries: parity, charge conjugation, and the discrete cubic 
rotations. The plateaux of the ratios (4) are typically found for t y G [10,15]. 3 The 
form factors Ti^^{q 2 ) are then extracted by minimising the x 2 on the corresponding set of 
plateaux of (4). When both mesons are at rest only the form factor T 2 can be computed, 
whereas in other kinematical situations we obtain all 3 form factors. In the following we 
will focus on Ti and T 2 . 

4 T B ^ K * (0) and T B ^ K * (0) /T 5 ^(0) 
4.1 Extrapolating to B at q 2 = 

As we already mentioned, in our lattice study we can extract the form factors at q 2 = 0, 
for each combination of kq-k^ in all three of our datasets (see table 1). The form factors 
that we directly compute on the lattice cover a range of q 2, s straddling around zero, so 
that either one of the kinematical configurations (9) coincides with q 2 = 0, or we have to 
interpolate the form factors calculated in the vicinity of q 2 = to q 2 = 0. In the latter 
case the results are insensitive to the interpolation formula used. 4 

A smooth linear mass interpolation (extrapolation) is needed to reach the H — > K* 
(H — > p) form factor, where H is our fictitious heavy-light meson that is accessible from 
our lattice. This is done by fitting to 

T^ v {q 2 = 0) = T^ v {q 2 = 0) = T H ^ V (0) = a H + (3 H m 2 P , (10) 

where mp is the light pseudoscalar meson, while an and Ph are the fit parameters. 
T H ^ K *(0) (T h -*p(0)) is then obtained after choosing m P = mf ys - (m* s '), where m x (vr) ph y s - 
are identified on the lattice by using the method of physical lattice planes [22]. Such 

3 Even after applying the available symmetries to the problem in hands, for each combination of HQ-n q 
we still have 73 correlation functions C^J^(t x ,t y ) when running over the ensemble of momenta (9). That 
means inspecting 4453 ratios (4) and from the corresponding plateaux we extracted 671 values for Ti j3 (<7 2 ), 
and 732 values of T 2 (q 2 ) form factor. We decided not to insert such formidable tables of numbers in this 
paper. A reader interested in those numbers can obtain them upon request from the authors. 

4 To check the insensitivity to the interpolation formula we used the forms discussed in eqs. (18,19) 
of the present paper, in addition to the pole/dipolc form, i.e., ?i(g 2 ) = T(0)/(1 — q 2 /inf) 2 , T 2 (q 2 ) = 
T(0)/(l-q 2 /m 2 ). 
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obtained values for T(0), together with 
ble 2. We also list our results for T H ~ 

(1 + m p ^ s -(3 H /a H )/(l + m^ s -(3 H /a H ). 



the masses in physical units, are given in ta- 
K * (0)/T H ^ p (0), which are simply obtained as 





«0 


M Hs [GeV] 


T H ^ K *(0) 




r—w—\ JJ , If* / „ \ 

T H K (0) 

T^^fOl 


QT^rp 1 


U. LZO 


1 70 (K\ 
1. < y i^oj 


n 7/1 
U. (Qyo ) 


U. ( U^o ) 


l.UO^O J 




u. 1ZZ 




U. ( y)(o) 




i.U 1 ^oj 




0.119 


2.29(6) 


0.65(6) 


0.60(8) 


1.09(4) 




U. 1 10 


9 ^Q^7"\ 




u.oo^y ) 


1. lU^O J 


Ofij 1 / 


U. l/o 


1.0 f 


U.oU^ll ) 


U. ( O^IOJ 


l.UO^O ) 




u. rzo 


l.o / \p) 


U.I 1 \o) 


n 7 r ^(^ o\ 
u. 1 oyLz j 








1 7"\ 
Z.Lo\t ) 


n 79 (7\ 


U.Oo^lU ) 


l.UO^O ) 




n 1 1 


O 10l(7\ 
Z.Oi) { 1 ) 


U.O ( (_ i ) 


U.oo^lU ) 


1 0^7^ 
l.UO^ 1 ) 




0.116 


2.62(8) 


0.62(8) 


0.57(11) 


1.09(9) 


Set 3 


0.1285 


1.80(6) 


0.75(7) 


0.72(10) 


1.02(3) 




0.125 


2.26(7) 


0.65(6) 


0.61(9) 


1.04(3) 




0.122 


2.62(9) 


0.57(5) 


0.51(7) 


1.07(4) 




0.119 


2.96(10) 


0.50(5) 


0.43(7) 


1.09(5) 




0.116 


3.28(11) 


0.44(5) 


0.38(7) 


1.11(6) 




0.114 


3.48(12) 


0.42(4) 


0.35(7) 


1.13(8) 



Table 2: The form factors at q 2 = computed directly on the lattice at a fixed value of the 
heavy quark and for all of our datasets. 

To extrapolate in the heavy quark mass we then use the heavy quark scaling law which 
tells us that T H ~* V (0) x mg 2 should scale as a constant, up to corrections proportional 
to l/mjj. Instead of the heavy quark mass, we may take the mass of the corresponding 
heavy-light meson consisting of a heavy Q-quark and the light s-quark. The reason for 
using the strange light quark is that it is directly accessible on our lattices whereas for 
the light -u/d-quark one needs to make an extrapolation which increases the error on the 
heavy-light meson mass. In other words we fit our data to 



T H - V (0) x MT S = Co + + 



r 3/2 



Cl 



c 2 



M. 



2 ' 



(11) 



where 00,1,2 are the fit parameters. From the plot in fig. 1 we see a pronounced linear 
behavior in 1/Mh s , which is why we will take the result of the linear extrapolation (02 = 0) 
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4.5 - 



> 

o 



3.5 - 



* 2.5 - 1_ 



X 



7.5 



□ H^Ky 

□ 5^^*y (quad.) 
■ B^Ky (lin.) 





0.2 



0.3 

7/M 



0.4 

[GeV 



0.5 



0.6 




Figure 1: In tie upper plot we show the extrapolation of the H — ► f£T*7 form factor (multiplied 
by Afg ; from 1/Mh„ directfy accessible on the lattice at (3 = 6.45, to 1/Mb s - Linear and 
quadratic fit to the data are denoted by the full and dotted lines respectively. The result of the 
quadratic extrapolation (empty square) is slightly shifted to left to make it discernible from the 
linear extrapolation result (filled square). The equivalent situation for the SU(3) breaking effect 
is shown in the lower plot. 
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as our main result. As it could be guessed from fig. 1, at f3 — 6.45, the extrapolated value 
does not change if we leave out from the fit the point corresponding to the lightest of our 
heavy quarks. Since we have more (and heavier) masses at (5 — 6.45, we prefer to quote 
the results obtained from that dataset (Set 3), namely, 

Ti£T**(0) = 0.24(3) , T£f (0) = 0.23(4) , (12) 

where "lin." and "quad." stand for the linearly and quadratically extrapolated form factors 
to 1/Mb s . The results of the strategy discussed in this subsection for all our lattices are 
listed in table 3. 









T£r K *(o)/T* n -;»(o) 


T q tT(0)/O(0) 


Set 1 


0.25(3) 


0.28(7) 


1.14(11) 


1.24(22) 


Set 2 


0.28(6) 


0.29(9) 


1.08(12) 


1.17(31) 


Set 3 


0.24(3) 


0.23(4) 


1.15(7) 


1.25(14) 



Table 3: The form factors resulting from the linear and quadratic extrapolation (11) as obtained 
from all of our 3 datasets specified in table 1. 



As an illustration, the linear fit with our data at (3 — 6.45 gives 

0.9(1) GeV 

72 * x " 



* (o) = 3 ' 6 ' 7 ' zf' 2 x 

MT 



(13) 



The slope in l/M Bs is very close to what has been observed in the lattice studies of the 
heavy-light decay constants [21], and of the B — > n semileptonic decay form factor (see 
eq. (19) in ref. [9]). 

From table 3 we see that the ratio of B — > K* and B — > p form factors has a large 
error. This error comes from T B ~* P (0), and in particular from the light mass extrapolation 
of the form factors to reach T H ^ P (0). That error is larger for larger m#, which gets further 
inflated after extrapolating to B. In contrast, the extrapolation to reach T H ^ K (0) is not 
needed as the K* mass falls in the range of the vector meson masses that are directly 
simulated on our lattices. When extrapolated to 5-meson, the SU(3) breaking ratio of the 
form factor has a large error and is very sensitive to the inclusion of the quadratic term in 
the extrapolation. 

Before closing this subsection, let us also mention that, contrary to HQET, one cannot 
match the short distance behaviour of our QCD results to the soft collinear effective theory 
(SCET) in which the T H ~* V (0) x mg 2 scaling law is manifest. This is still an unsolved 
theoretical problem and hopefully a recent development based on ref. [23] will help solving 
it. We note however that the inclusion of the matching of the tensor current anomalous 
dimension of QCD with HQET produces a numerically marginal effect (1 -j- 2% on the 
central values). We hope that the similar will hold once such a matching of QCD with 
SCET becomes possible. 



9 



4.2 Extrapolating to B at q 2 ^ and then to q 2 = 

To check on the results obtained in the previous subsection we now also employ the standard 
method and extrapolate our results at fixed v-p' to Mb- The main assumption here is that 
the HQET scaling laws are valid for all our v-p', i.e., not only for those that are very small 
compared to the heavy meson mass. 

• From our directly accessed masses and g 2 's, one identifies vp' = (Mfj+rriy — q 2 ) / '2M h , 
where H is the heavy-light meson and V stands for either K* or p. In physical units, 
the range of available v-p' is nearly equal for all of our lattices, namely 0.9 GeV < 
v-p' < 1.8 GeV. We emphasize that the kinematical configurations which we were 
able to explore are those listed in eq. (9). Proceeding like in ref. [9], we chose 5, 
6, 7 equidistant v-p' points for our dataset 1, 2, 3, respectively. The form factors, 
T\ t2 (v-p'), are then linearly interpolated (extrapolated) to rriK* (jn P ) for each of our 
heavy quarks. 

• We construct 

$ 1 (M H ,v-p')=w(M H )^=fi, $ 2 (M H ,vp')=w(M H )T 2 (vp') x /M^, (14) 

V M H 

which, in HQET, are expected to scale as constants, up to corrections oc 1/Mj^. 
Fitting our data to 

*,, 2 (MW)=^) + ^ + ^, (15) 

either linearly (d 2 = 0) or quadratically (d 2 ^ 0), we can extrapolate to the £?-meson 
mass. Since we use the scaling law which is manifest in HQET the factor w(M H ) 
in eq. (14) accounts for the mismatch of the leading order anomalous dimensions in 
QCD ('Jt = 8/3) and in HQET (7V = —4) for the tensor density, namely, 

w{Mh) = ) x Uw*) ) 

The numerator in the first factor match our QCD form factors with their HQET 
counterparts, while the dominator does the opposite to the result of the extrapolation 
to M B . The second factor, instead, provides the NLO evolution from /i — 1/a to 
H = M H . For N f = 0, fa = 11, and J T = 2.53. 

• The extrapolation (15) is made both linearly and quadratically. The differences 
between the corresponding results are essentially indistinguishable, as it can be seen 
in fig. 2 where we plot the results for B — > K*^* form factors (7* stands for the 
off-shell photon) obtained at /3 — 6.45 (Set 3). These and the similar results we 
obtained at /3 — 6.2 are collected in table 4. 

To reach the physically interesting case of the photon on-shell one needs to assume some 
functional dependence of the form factors and extrapolate the results of table 4 down to 



J, 



a s {l/a) - a s {M H ) 

47T 



• (16) 
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Figure 2: The form factors Ti^(q 2 ) relevant for B — > K*~/* decay, obtained after extrapolating 
(linearly and quadratically) our data at (5 = 6.45 in inverse heavy meson mass. Also shown are 
the curves fitting the q 2 dependence to the expressions given in eqs. (18,19). 
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Set 1 



Set 2 



q 2 [GeV 2 ] 


Tr K *(q 2 ) 




12.2 
13.6 
15.1 
16.5 
17.9 


0.76(24)1;°° 
0.94(22)1;°° 
1.19(19)1;°° 
1.53(17)1;°° 
2.03(21)1;°° 


0.50(13)1;°° 
0.54(12)1;°° 
0.60(10)1;°° 
0.64(6)1;°° 
0.72(2)1;°° 









q 2 [GeV 2 ] 




T^ K \q 2 ) 


13.5 
14.4 
15.3 
16.1 
17.0 
17.8 


0.70(17)l;g 

o.8i(i8)i;55 

0.95(19)l;g 
1.11(21)1;° 3 
1.32(26)1;°° 
1.56(36)1;°° 


0.52(9)1;°° 
0.56(9)1;°° 
0.61(9)1;°° 
0.66(9)1;°° 
0.72(9)1;°° 
0.80(9)1;°° 



Set 3 



q 2 [GeV 2 ] 


Tr K \q 2 ) 


T 2 B - K *{q 2 ) 


11.2 
12.4 
13.6 
14.7 
16.0 
17.2 
18.3 


0.73(10)1;°° 
0.86(11)1;°° 
1.02(13)1;°° 
1.23(16)1;°° 
1.52(23)1;2° 
1.93(34)1;°° 
2.54(55)1;°° 


0.39(6)1;°° 
0.43(6)1;°° 
0.47(6)1;22 
0.52(6)1;2° 

o.58(7)l;22 

0.66(8)1;2° 
0.76(10)1;°° 



Table 4: The values of the B — ► K*j* form factors at several values of q 2 . The first error in each 
result is statistical and the second is the difference between the results of linear and quadratic 
extrapolation in 1/Mh to I/Mb- 
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q 2 = 0. It is very easy to convince oneself that the form factors 7\ and T 2 satisfy the 
constraints very similar to those that govern the shapes of F + and F semileptonic heavy 
to light pseudoscalar form factors. More specifically: 

o The nearest pole in the crossed channel, J p = 1~ , which contributes to the form 
factor Ti, is Mb* = 5.42 GeV in B — > K* transition. The form factor T 2 , instead, 
receives the contribution from heavy J p = 1 + resonances and multiparticle states 
both below and above the cut (M B + my) 2 . 

o HQET, which is relevant to the region of large g 2 's, suggests that the form factors 
scale with heavy quark/meson mass as Ti ~ \[M and T 2 ~ l/VM [11], and therefore 
both form factors cannot be fit to the pole-like shapes. 

o For the high energy region of the light meson in the rest frame of the heavy (q 2 — > 0), 
we also have the scaling laws Ti(E) ~ \[M / E 2 ~ M~ 3 / 2 . Similar holds for the T 2 (E) 
form factor, i.e., both form factors scale as M~ 3 / 2 [16,17]. Moreover, the two are 
related via 

Ti(E) — ^-T 2 (E) . (17) 



Thus the situation is analogous to the one in B — > irlv decay, and we can use the simple 
parameterisation of ref. [24], 

T (a 2 ) - T{0) T (a 2 ) - T(0) (IR) 

Tl{q) - (l-^)(l-a^)' Uq ) -l^WW' { ) 

where q 2 = q 2 /M|». If we relax the Xi/T 2 = M/2E constraint, then a simple form (18) 
becomes 



Our data from table 4 cannot distinguish between the two of the above parameterisations 
and in both cases we end up with the same value for T(0). As in the previous subsection, 
as our final results we will quote those obtained at (3 = 6.45, for which more and heavier 
mesons were accessible, 

T£T K \0) = 0.22(3) , T*Z K \0) = 0-24(4) . (20) 
The results for all three datasets are collected in table 5. 



5 Final results and conclusion 

As our final results we will quote those obtained at (3 — 6.45, since they have smaller 
discretisation errors. As central value we take the results obtained from the first method 
(subsec. 4.1) which are given in table 3 at \i = 1/a, in the (Landau) RI/MOM scheme. 
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rpB-+K* (f\\ 
J lin. l U J 


rpB-^K* (C\\ 
J quad. l U J 


J lin. l U J/ J lin. V 3 ) 


J quad. l U J/ J q U adA U J 


Set 1 


0.25(6) 


0.28(10) 


1.2(3) 


1.3(7) 


Set 2 


0.23(6) 


0.23(6) 


1.1(3) 


1.1(2) 


Set 3 


0.22(3) 


0.24(4) 


1.3(2) 


1.3(3) 



Table 5: The results of the extrapolation of the form factors from table 4 assuming the q 2 - 
dependence given in eqs. (18,19). 



The matching to the MS scheme is 1 at 2-loop accuracy, but we still need to run our result 
from fi = 1/ a to nib which is made via 



T B ^ K *(0;m b )=T(0;l/a) 



a s (l 



-7t/2/3o 



a s (m b ) J 



I- J, 



a s (l/a) - a s {m b ) 

47T 



(21) 



For Nf = and m b = 4.6(1) GeV, the running factor is very close to 1, i.e., 0.99 (0.97) for 
the form factors computed at (5 — 6.45 (f3 = 6.2). Our final result is 



T B^K* (g 2 = Q. ^ = ^ = Q2A ± Q.03+0.04 



(22) 



The spread of central values presented in table 3 at (5 — 6.2 (when multiplied by 0.97) 
are used to attribute a systematic error to our final result (22). The values in table 5 
are already at the scale \x = M B m b) and those given in eq. (20) are fully consistent 
with (22). Our number is smaller when compared to the QCD sum rule results, 



T (0) = 0.33(3) [7], 0.38(6) [8], 



(23) 



although the recent progress show the tendency of lowering the ccentral value obtained by 
using LCSR, i.e., T B ^ K *(0) = 0.31(4) [25]. Concerning the previous lattice studies [26], 
most of them were made at the time before the T(0)M^ //2 -scaling law was known or before 
a significant statistical quality of the lattice data was feasible. 

Our values for the ratio T B ~> K * (0)/T B ^f>(0), on the other hand, are unstable, and we 
will only quote the average of the values obtained by using two methods discussed in this 
paper at f3 — 6.45, 



T B ^*(0)/T^(0) = 1.2(1), 



(24) 



as our best estimate. This result agrees with the most recent LCSR estimate, T B ~* K * (0)/T B ^ f 
1.17(9) [27]. 

The method employed in this work relies on the use of a propagating heavy quark. In 
order to reach smaller values of g 2 's in an effective theory of heavy quark the so called 
moving NRQCD has been developed, but the numerical quality of the signal does not 
appear very encouraging so far [28]. Clearly, our method cannot be used to obtain a very 
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accurate value of T(0), mainly because the heavy quark extrapolations are involved. If 
we are to work with the physical 5-quark mass on the lattice, we would need a very small 
lattice spacing, i.e., at least a" 1 ~ 10 GeV. For the volume corresponding to the lattice 
size La = 2 fm that would require the simulations on the lattice with 100 3 spatial points. 
Moreover, on the lattice with the periodic boundary conditions, the q 2 = point is reached 
when the energy of the vector meson (in the rest frame of the heavy) is 

\H\ = f K ' xLa= (2.07 fm" 1 ) x La . (25) 

Therefore on a lattice of the size La = 2 fm, one needs to give the kaon a momentum \n\ ph 4, 
for which it is highly unlikely to observe any signal of the correlation functions (5). Even if 
one organises the kinematics so that the momenta are shared between B- and i^*-mesons, 
the required spatial momenta are still too large for the reasonably accurate computation 
of the correlation functions (5). Therefore, a progress in improving the quality of signals 
of correlation functions when the spatial momenta \n\ > 1 would be highly welcome. 
Summarising, as of now, a significant improvement on the precision of the B — > K* form 
factors does not look promising even in quenched approximation. The methodology of 
the extraction of the form factors can, however, be improved by combining the correlation 
functions (5) in double ratios, in a way similar to what has recently been implemented in 
the lattice computation of heavy— dieavy [29] and light— dight pseudoscalar meson decay 
form factors [30] . 
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